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Abstract. Accreting X-ray pulsars are among the most luminous objects in the X-ray sky. In 
highly magnetized neutron stars (B ~ 10 12 G), the flow of matter is dominated by the strong 
magnetic field. The general properties of accreting X-ray binaries are presented, focusing 
on the spectral characteristics of the systems. The use of cyclotron lines as a tool to directly 
measure a neutron star's magnetic field and to test the theory of accretion are discussed. We 
conclude with the current and future prospects for accreting X-ray binary studies. 
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1. Introduction 

Accretin g X-ray binaries were discovered in 
the 70 's dGiacconi et al.ll97ll lTananbau m et alj 
Il972h and are among the most luminous ob- 
jects in the X-ray sky. Forty years after their 
discovery, more than 400 X-ray binaries have 
been observed, and there is today a deeper 
knowledge of the physical processes involved 
in the accretion of matter onto the compact 
object. In X-ray binaries, matter is accreted 
from a donor star onto a compact object (white 
dwarf, neutron star or black hole). X-ray emis- 
sion originates as a result of the conversion of 
the gravitational energy of the accreted mat- 
ter into kinetic energy. According to the na- 
ture of the donor star, X-ray binaries can be 
classified as High-Mass X-ray Binaries (here- 
after HMXBs) or Low-Mass X-Ray Binaries 
(hereafter LMXBs). Typically, HMXBs have 
young optical companions of spectral type O 



or B and mass M > 5M , and high mag- 
netic fields B ~ 10 12 G. LMXB systems have 
older optical companions, with masses in gen- 
eral M < 1 M Q , and lower magnetic fields B ~ 
jq9-10q X-ray binaries are numerous objects 
in the Galaxy, with 1 14 HMXBs in the Galaxy 
and 128 HMXBs in the Magellanic Clouds, 
and 187 LMXBs in the Galaxy and M agellanic 
Clouds dLiu et al J2005Ll2006Ll2007h . In accret- 
ing neutron stars, due to the inclination of the 
neutron star's rotational axis with respect to the 
magnetic axis, the X-ray emission appears pul- 
sated to a distant observer. 

In this review, we concentrate on the prop- 
erties of X-ray binaries with a highly magne- 
tized neutron star (B ~ 10 12 G). The basic as- 
pects of accreting X-ray pulsars are described, 
and the current status of the observations and 
theory of these systems is presented. 
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2. Accretion onto the neutron star 

The accretion of matter in X-ray binaries can 
take place via Roche lobe overflow or via direct 
wind accretion. Roche Lobe overflow takes 
place when the donor star evolves, expands and 
fills its Roche lobe. The matter exceeding the 
Roche lobe is no longer gravitationally bound 
to the star and can be accreted by the compact 
object through the inner Lagrange point. The 
accreted matter has a large amount of angular 
momentum, and therefore forms an accretion 
disk around the compact object instead of be- 
ing accreted directly. This typically occurs in 
LMXBs. 

Wind accretion takes place in binary sys- 
tems with a O or B donor star, namely HMXBs. 
These stars have very strong stellar winds, with 
a mass loss rate that can be as high as M ~ 
10~ 4 - 10~ 6 M o /yr. The compact object typi- 
cally orbits the donor star at a distance of less 
than one stellar radius, and therefore the com- 
pact object is deepl y embeded in the stellar 
wind (e.g., Vela X-l. lNagase etaDl986h . 

A large fraction of X-ray binaries (about 
60% of the HMXBs) falls in the category of 
Be/X-ray binary systems. These systems have 
quite eccentric orbits, and their X-ray activity 
is related to the presence of a disk of material 
that surrounds the equator of the Be star. They 
can have giant (type II) or normal (type I ) out- 
bursts, and long quiescence periods. See iReigl 
d201 lh for a recent review of Be/X-ray bina- 
ries. 

Independently of the way matter is trans- 
ferred from the donor star to the compact ob- 
ject, at a distance close to the compact ob- 
ject the flow of matter is dominated by the 
strong magnetic field of B ~ 10 12 G. The in- 
falling matter follows the magnetic field lines, 
forming accretion columns on the neutron star 
magnetic poles. The accreting plasma couples 
to the magnetic field lines at the Alfven ra- 
dius r mag = 2.9 x 10 8 Mj /7 fl 6 2/7 L3 7 2/ y /7 cm, 
where M\ is given in units of 1 M Q , R(, is 
given in units of 10 6 cm and L37 is the luminos- 
ity in units of 10 37 ergs -1 . For typical neutron 
star parameters, r mag ~ 1800 km. Accretion 
can take place as long as the Alfven radius 
is smaller than the co-rotation radius r co = 
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Fig. 1. Schematic view of the accretion col- 
umn. Seed photons are created along the column 
via bremstrahlung and cyclotron emission, together 
with black body seed photons emitted from the sur- 
face of the thermal mound. (Figure 1 from Becker 
& Wolff d2007h " reproduced by permission of the 
AAS). 



(^jr) ' (defined as the radius at which the an- 
gular velocity of the magnetosphere u>r and 
the Keplerian velocity of the disk y/GM/r are 
equal). If this is not the case, the propeller 
effect can inhibit the accretion (Illarionov & 
Sunvaev Il975h . Such an effect could for in- 
stan ce explain the off-states ob served in Vela 
X-l dKrevkenbohm et alJl2008h . 

The luminosity produced by an accreting 
object is given by L = GMm/R, where M, 
R are the mass and radius of the compact 
object, and m is the mass accretion rate. As 
will be further discussed in Sec. |4] there is a 
limit on the luminosity that an accreting body 
can have, the Eddington luminosity, obtained 
by balancing gravitational force and radiation 
pressure. For a spherically symmetric accretion 
and a steady flow, the Eddington luminosity is 

LEdd= *^ wL3xl0 38 M ergs -l i 

3. Formation of the spectral 
continuum 

The X-ray spectrum is typically well described 
by a power law with an exponential cutoff. 
Numerous authors have made attempts to de- 
rive the shape of X-ray pulsa r spectra an- 
alytically or numerically (e.g., iNagell Il98ll 
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Meszaros & Nagellll985llBurnard et a"Dll99lL 
Becker & Wolffl bo05h . However no self- 
consistent, general model applicable to X-ray 
sources has been established, due to the com- 
plexity of the physical processes that take place 
in the accretion column and in the magneto- 
sphere. The radiatio n spectrum or "standard" 
X-ray continuum dWhite et al.ll983h is a power 
law in the ~ (5 — 20) keV energy range with an 
exponential cutoff at energies ~ (20 - 30) keV 
(ICoburn et al.l [2002). Additionally, iron fluo- 
rescence Fe Ka lines are produced in the cir- 
cumstellar material, and cyclotron absorption 
lines also modify the continuum emission (see 
Sec. |D. 
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Fig. 2. Theoretical spectrum of Her X- 1 . The to- 
tal spectrum (red line) and individual components 
are shown, together with a BeppoSAX spectrum of 
the source (black circles). (Figure 6 from Becker 
& Wolff d2007h .' reproduced by permission of the 
AAS). 
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Fig. 3. Unfolded spectrum of 4U 01 15+64 and best 
fit model. The solid line shows the best fit model, 
sum of the dashed line that represents the thermal 
Comptonization, the dot-dashed line that shows the 
column emission, and the dot-dot dashed line that 
shows a narrow Gaussian line. Low er panel: resid- 
uals fr om the best fit model. From IFerrigno et"aTI 
(2009), reproduced with permission © ESO. 

The mod el is available for spe ctral fitting 
in XSPEC (IFerrigno et al.l 120091) . The au- 
thors have applied the model to BeppoSAX 
observations of the accreting Be/X-ray binary 
4U 01 15+64 during an outburst. From this 
study, Ferrigno et al. found that the emission 
above ~1 keV could be due to thermal and bulk 
Comptonization of the seed photons produced 
by cyclotron cooling of the accretion column, 
plus emission at lower energies due to ther- 
mal Comptonization of a blackbody compo- 
nent from a diffuse halo close to the neutron 
star surface. The data, best fit model, and resid- 
uals are shown in Fig. [3] 



Recently. iBecker & Wolj (120071) have de- 
veloped a model for the continuum formation 
to reproduce the spectra of bright accreting X- 
ray pulsars, for which a radiative shock located 
above the stellar surface dominates the spectral 
formation. The total observed spectrum is a re- 
sult of the bulk and thermal Comptonization 
of bremstrahlung, black body and cyclotron 
seed photons. A schematic view of the model 
is shown in Fig. Q] In Fig. [2] a comparison 
between the theoretical spectrum of Her X-1 
and a BeppoSAX observation of the source is 
shown. 



4. Cyclotron Resonance Scattering 
Features 

The energy of the electrons in the highly mag- 
netized plasma is quantized into Landau lev- 
els. Due to the resonant scattering of photons 
off electrons, absorption-like features can be 
formed in the X-ray spectrum of X-ray pul- 
sars (cyclotron resonance scattering features or 
simply cyclotron lines). The energy spacing 
between the Landau levels is given by 

eB 

E cyc =/z = ll.6keV-.Bi2 (1) 



Caballero, I., Wilms, J.: X-ray pulsars: a review 



285 



Tabl e 1. List of sources with cycl otron line(s) s ignifi cantly detected in their spectrum, updated 
from Ma kishima et al.l (11999 ') and He indl et al.l (120041) . The cyclotron line energy, the discovery 
instrument and reference are given. For the last five sources listed there is weak evidence for a 
cyclotron line but this has not been firmly confirmed. 



Source 


£n,cyc (keV) 


Discovery 


Swift J 1626.6-5 156 


10 


RXTE 


4U 01 15+64 


14,24,36,48,62 


HE AO- 1 


4U 1907+09 


18,38 


Ginga 


4U 1538-52 


20 


Ginga 


VelaX-1 


24,52 


HEXE 


V 0332+53 


26,49,74 


Ginga 


CepX-4 


28 


Ginga 


CenX-3 


28.5 


RXTE 






BeppoSAX 


XPer 


29 


RXTE 


MXB 0656-072 


36 


RXTE 


XTEJ1646+274 


36 


RXTE 


4U 1626-67 


37 


RXTE 






BeppoSAX 


GX 301-2 


37 


Ginga 


HerX-1 


41 


Balloon 


A 0535+26 


46,100 


HEXE 


1A 1118-61 


55 


RXTE 


EXO 2030+375 


11? 


RXTE 




63? 


INTEGRAL 


GS 1843+00 


20? 


Ginga 


OAO 1657-415 


36? 


BeppoSAX 


GROJ1008-57 


88? 


CGRO 


LMCX-4 


100? 


BeppoSAX 



Reference 



DeCesar et al. (2009) 
Wheaton etal. (19791 
Makis hima & Mihara (1992) 
Clark etal. (1990) 
Kendziorra et al. (1992) 
Makishima et al. ( 1 990) 
Miha ra et al. (1991) 
Heindl & Ch akrabartv (1999) 
Santangelo et al. (1998) 
Coburn et al. (2001) 
Heindl et al. (2003) 
Hein dl et al. (2001) 
Heindl & Ch akrabartv (1999) 
Orlandini et al. (1998) 
Makishima & Mihara (1992) 
Triim per et al. (1977) 
Kendziorra et al. ( 1992) 
Doroshenko et al. (2010) 
Wilson et al. (20081 
Klochkov et al. (2008a) 
Mihara (1995) 
Orlan dini et al. (1999) 
Shraderetal. (1999) 
La Barbera et al. (2001) 



where B12 is the magnetic field in units of 10 12 
G. Therefore, the measurement of a cyclotron 
line in the spectrum of a highly magnetized ac- 
creting neutron star provides a direct measure- 
ment of the neutron star's magnetic field. Due 
to the strong gravitational field around the neu- 
tron star, the cyclotron line energy is gravita- 
tionally red-shifted, and the observed line en- 
ergies have to be corrected as 

E n = nE cyc = (1 + z)£„ )0bs (2) 

where the gravitational redshift is z ~ 
1 .25 ... 1 .4 for typical neutron star parameters. 

Since the discovery of a cyclotron line 
in the X-ray spectrum of Her X-l in 1976 
dTruemper et al.l fl978h . cyclotron lines have 
been observed in 16 accreting pulsars (see 
e.g. [Makishima et all I1999L iHeindl et alJ 
2004, for reviews, and Table [T] for a list of 
the currently known sources that exhibit cy- 
clotron lines). Several sources exhibit a fun- 
damental line plus harmonics, like V 0332+53 



dKrevkenbohm et aT]|2005l) . The record holder 
is 4U 01 15+64, show ing up to five cyclotron 
lines in its spectru m dSantangelo et al.lll999l 
IHeindl et al.lfl999l) . 

Cyclotron lines are extremely powerful 
tools since they provide the only direct way to 
determine the magnetic field of a neutron star, 
and they can also probe the change of the ac- 
cretion structure and plasma properties with lu- 
minosity. 

Matter is channeled onto the magnetic 
poles by the magnetic field lines creating ac- 
cretion funnels, where the X-ray emission is 
ori ginated. As proposed in the pioneering work 
of Basko & Sunvaev (119761) . there are two 
different accretion regimes according to the 
mass accretion rate (and hence the luminos- 
ity). For luminosities higher than a critical lu- 
minosity, which is general ly comparable to th e 
Eddington luminosity Le dNelson et all 1993). 
a radiation dominated shock is formed in the 
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Fig. 4. Evolution of the cyclotron line energy 
of V 0332+53 with th e X-ray luminosity. From 
iTsveankov et al.l d2010h . reproduced with the per- 
mission from MNRAS. 



accretion column. In this case a "fan beam" 
emission pattern is predicted, in which pho- 
tons escape perpendicular to the accretion col- 
umn. For luminosities lower than the critical 
luminosity, matter can fall almost freely until 
the neutron star surface, causing the braking of 
the plasma by a hydrodynamical shock close 
to the neutron star surface. In this case, X-rays 
are able to escape vertically along the accre- 
tion column, giving rise to a "pencil beam" pat- 
tern. For high luminosity sources it is gener- 
ally believed that an increase in the mass ac- 
cretion rate (or increase in luminosity) causes 
an inc rease in the height of the accretion col- 
umn dBurnard et al.lll99ll based on Basko 
& Sunyaev 119761) . Assuming a dipole mag- 
netic field, the height of the accretion column 
can be estimated in terms of the cyclotron en- 
ergy. Therefore, cyclotron line energy mea- 
surements allow us to test this theory. This 
relation was experimentally confirme d with 



Ginga observations of 4U 1 15+6 4 (Mihara 



1995L iMihara et all I2004L [Nakajima et al 



2006) and RXTE observations of 4U 01 15+64 
dNakaiima et alJl2006llTsv gankov et alj|2010l) 



and V 0332+53 (Tsvgank ovet al.ll2007h . For 
these sources, a negative correlation between 
the X-ray luminosity and the cyclotron line en- 
ergy has been observed (see Fig. |4). 

For low luminosity sources in the sub- 
Eddington regime, a different behavior is ex- 




max. ASM cls/s 



Fig. 5. Dependence of cyclotron line energy with 
the X-rav luminositv for HerX-1. From Staubert 
et al. (2007), reproduced with permission © ESO. 



pected. Staube rtet al.l d2007h have discovered 
for HerX-1 a positive correlation between the 
cyclotron line energy and the luminosity (see 
Fig. [5]). They find a change in the cyclotron line 
of 5% for a luminosity variation of AL/L ~ 
1. In this case, no radiation dominated shock 
is expect ed to form in the accret ion column. 
Based on B asko & Sunvaevld 19761) and Nelson 
et al. (1993), where the physics of accretion 
in low-luminosity sources is studied, Staubert 
et al. d2007l) find the following relation be- 
tween the cyclotron line energy and the lumi- 
nosity: AE cyc /E cyc — 3^^ where Z, is the 
height of the scattering region above the neu- 
tron star surface and R is the neutron star ra- 
dius. For AL/L ~ 1, this implies a change of 
~3%, very close to the 5% obtained from the 
HerX-1 observations. 

Contrary to the above sources, mea- 
surements of the cyclotron line energy of 
A 0535+26 during a normal outburst in 2005 
show no significant correlation between the cy- 
clotron line energy and the luminosity. This 
suggests that the line forming region does not 
chan ge with the mass accr etion rate of the sys- 
tem dCaballero et al.ll2007l) . A question that is 
still under debate is whether the transition be- 
tween the two regimes has been observed. 

Changes in the cyclotron line energy have 
been reported not only with the luminosity, but 
also with the pulse phase. The cyclotron line 
energy values discussed above correspond to 
phase averaged spectra. However, due to the 



Caballero, I., Wilms, J.: X-ray pulsars: a review 



287 



30 
12 
10 

8 

f 6 
4 

2 

1.4 
1.2 
1.0 
0.8 
0.6 
0.4 



0.0 0.2 0.4 0.6 0.8 
Pulse phase 



0.0 0.2 0.4 0.6 0.8 1.0 
Pulse phase 



Fig. 6. Variation of the cyclotron line parameters 
with the pulse phase for Her X-l, during the start of 
the main-on state (Jeft ) and at the end of the m ain- 
on state (right). From Kloc hkov et alJ J2008bl) . re- 
produced with permission © ESO. 

rotation of the neutron star, during a rotational 
phase the observer is looking into different re- 
gions of the neutron star surface and/or ac- 
cretion column, and therefore different mag- 
netic fields and continuu m parameters can be 
expect ed. For GX 301-2, Krevkenboh m et ail 
(2004) have shown that the cyclotron line en- 
ergy varies by 25%, implying a variation of the 
magnetic field betwe en 3.4 X 10 12 and 4.2 X 
10 12 G. For Her X-l, iKlochkov et ail d2008bl) 
have also found a variation of the cyclotron 
line and continuum parameters with the pulse 
phase (see Fig. [6j>, that they interpret with dif- 
ferent viewing angles of the neutron star along 
the rotational phase. 

Until recently, only phenomenological 
models (Gaussian or Lorentzian curves) were 
used to describe the cyclotron lines observed 
in the spectra of accreting X-ray pulsars. 
However, Monte Carlo simulations of the prop- 
agation of photons through a low density 
plasma assumed to be threaded by an uniform 
magnetic field have revealed that the cyclotron 
lines are expected to vary in shape, depth and 
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Fig. 7. Simulated cyclotron line shapes for 
B/Ban = 0.04 as a function of angle (decreasing 
upwards) and optical depth. A cylindrical geometry, 
a plasma temperature of kT = 3 keV, Thomson opti- 
cal depths of 3 x 10 4 (left) and 3 x 10 3 (right) are as- 
sumed. Line shapes are rather complex and strongly 
depen d on the physica l and g eometrical parameters. 
From lSchonherr et al.l (2007). reproduced with per- 
mission © ESO. 



width over the pulse (lArava & Hardindll999[ 
lArava-Gochez & Harding! l2000l Schonherr 
et al.l2007, see Fig. [7]). A complex line pro- 
file has also em erged from observa t ions ( e.g., 
in V 0332+53, iPottschmidt et alJ l2005n . A 
Monte Carlo code based on the code of Araya 
&Harding d 19991) . lArava-Gochez & "H arding 
(2000) has been impl emented in XSPEC by 
dSchonherr et al.l2007l) . Fi rst fits have been per - 
formed to Cen X-3 data (ISuchv et al.ll2008h . 
The new model allows to narrowly constrain 
the main physical parameters of the system, as 
the plasma geometry (slab versus cylinder), the 
electron parallel temperature, and density. The 
latest version of the model allows for complex 
geometries, ma gnetic field grad ients, and ve- 
locity gradients (lSchwarmll2010h . 



5. Conclusions 

We have reviewed the general properties of ac- 
creting X-ray pulsars, in particular the spec- 
tral formation and cyclotron lines. Many as- 
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pects of accreting X-ray pulsars have not been 
discussed in this review, like pulse profile for- 
mation. Modeling of pulse profiles has been 
bv jWang & Welter 



perfor med for instance by IWang & 



ry [Wans 

(1981). lMeszaros & NapeM (Il985l) . and lLeahvl 
dl99lh . An alternative pulse profile decom- 
position method has been proposed bv Kraus 
et al. dl995h . that allows to infer geometri- 
cal parameters of the neutron star, and to re- 
construct the emission beam pattern that can 
be compared to model predictions (see Kraus 
et al. (120031) and e.g. lSasaki et al.ll2010l)' . Also 
not discussed in this review are the relativisti- 
cally broadened lines that have been observed 
in several neutron star X-ray binaries. These 
can be used to constrain neutron star parame- 
ters and the equation of state, see, e.g., Cackett 
et al J2008I) but also iDone & Diaz Trigol 
d2010h . 

Forty years after the discovery of the first 
accreting neutron star, the continuum forma- 
tion is now better understood, but still many 
questions are open. There are now 16 neutron 
stars with direct magnetic field determinations 
- now enough data available for both, indi- 
vidual studies, and study as a class. The cy- 
clotron line behavior with the X-ray luminosity 
allows deeper probing of accretion column the- 
ory. Good numerical models for cyclotron line 
formation are available, and the behavior of the 
lines is roughly in agreement with predictions 

of Monte Carlo computations; 

In the future, M AXI dMatsuoka et alJ 



120091) and Swift/BAT dBarthelmv et al.ll2005h 
will continue to monitor the X-ray sky and 
discover neutro n star outbursts. ASTROSAT 
(lO'Brienl I20H will allow routine measure- 
ments of broadband spectral shape and mag- 
netic fields, NuStar dHarrison et al.ll2005l) will 
start to resol ve the shape of cyc lotron lines, 
while LOFT dFeroci et al] 1201 Oh and Athena 
will allow studies at characteristic variability 
timescale. With the current and future mis- 
sions, we have reached a very exciting era in 
which high quality data are available together 
with physical models that will allow a deeper 
understanding of the extreme physics involved 
in the accretion of matter onto a neutron star. 
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